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New critical behavior in unconventional superconductors and superfluids is established and de- 
scribed by the Wilson-Fisher renormalization-group method. For certain ordering symmetries a new 
type of fluctuation-driven first order phase transitions at finite and zero temperature are predicted. 
The results can be applied to a wide class of ferromagnetic superconducting and superfluid systems, 
in particular, to itinerant ferromagnets as UGe2 and URhGe. 
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I. INTRODUCTION 

In this paper new critical behavior in unconventional 
ferromagnetic superconductors and superfluids is estab- 
lished and described. This behavior corresponds to an 
isotropic ferromagnetic order in real systems but does not 
belong to any known universality class [l| and, hence, 
it could be of considerable experimental and theoreti- 
cal interest. Due to crystal and magnetic anisotropy, a 
new type of fluctuation-driven first order phase transi- 
tions occur, as is shown in the present investigation. The 
novel fluctuation effects can be observed near finite and 
zero temperature ("quantum") phase transitions 1, 2] in 
wide class of ferromagnetic systems with unconventional 
(spin-triplet) superconductivity or superfluidity. 

The present study has been performed on the special 
example of intermetallic compounds UGe2 and URhGe, 
where the remarkable phenomenon of coexistence of itin- 
erant ferromagnetism and unconventional spin-triplet su- 
perconductivity 3] has been observed M. For example, 
in UGe2, the coexistence phase occurs [4!] at tempera- 
tures < r < 1 K and pressures 1 < P < Pq ^ 
1.7 GPa. A fragment of {P,T) phase diagrams of itin- 
erant ferromagnetic compounds 4\ is sketched in Fig. 1, 
where the lines Tp{P) and Tc{P) of the paramagnetic(P)- 
to-ferromagnetic (F) and ferromagnetic-to-coexistence 
phase (C) transitions are very close to each other and 
intersect at very low temperature or terminate at the 
absolute zero {Po,0). At low temperature, where the 
phase transition lines are close enough to each other, 
the interaction between the real magnetization vector 
M{r) = {Mj{r);j = l,...,m} and the complex order 
[ parameter vector of the spin-triplet Cooper pairing 
' V(r) = {Mr) = 1, ••••n/2} (n = 6) can- 

not be neglected [l| and, as shown here, this interaction 
produces new fluctuation phenomena. 

Both thermal fluctuations at finite temperatures (T > 
0) and quantum fluctuations (correlations) near the P~ 
driven quantum phase transition at T = should be 
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FIG. I: {P,T) diagram with a multicritical point {Po,T 
0) at very low temperature. Para- (P), ferromagnetic (F), 
and coexistence (C) phases, separated by the lines Tf{P) and 
Tc{P) of P-F and F-C phase transitions, respectively. 



considered but at the first stage one may neglect the 
quantum effects [2] as irrelevant to finite temperature 
phase transitions {Tp ^ Tc > 0). The present treat- 
ment of recently derived free energy functional Q by the 
standard Wilson-Fisher renormalization group (RG) fl| 
shows that unconventional ferromagnetic superconduc- 
tors with isotropic magnetic order (m = 3) exhibit a 
quite particular multi-critical behavior for any T > 0, 
whereas the magnetic anisotropy (to — 1,2) generates 
fluctuation-driven first order transitions (see Sec. Ill and 
IV). 

As shown in Sec. II certain terms in the general free 
energy Hamiltonian [6| , which are relevant for the mean- 
field analysis @ become irrelevant within the RG 
framework. This leads to a considerable simplification 
of the RG treatment but on the other side specific sym- 
metry properties of the relevant Hamiltonian terms make 
the same RG treatment quite nontrivial (Sec. III). The 
RG study is performed by the known e-expansion (see, 
e.g., Ref. 0). 

A particular feature of the present theoretical consid- 
eration is the breakdown of usual e-expansion in the one- 
loop approximation for RG equations (e = 6 — d). Thus, 
the latter cannot yield conclusive results and one is forced 
to derive the two-loop RG equations. The higher order 
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of the theory restores the e-expansion in a modified form 
- expansion in non-integer powers (e^^^, e"'^/^, . . . ) of e. In 
the framework of two-loop RG equations one obtains the 
above mentioned new critical behavior. 

This theoretical scenario looks quite similar to the 
breakdown of the usual e-series and the e-expansions in 
non-integer powers of e = (4 — d), known for the first 
time from the theory of critical phenomena in anisotropic 
disordered systems (see, e.g., Ref. [9]). It should be 
emphasized that the mentioned similarity between the 
present theoretical analysis and that in certain disor- 
dered systems cannot be extended beyond some general 
features of the e-expansions. The mechanism leading to 
the e = (4 — d)-expansion in non-integer powers of e in 
anisotropic disordered systems is quite different from the 
mechanism revealed here. The latter is a result of the 
particular symmetry of the interaction between the fields 
ip and M. 

In order to re-derive the results in Sec. Ill, one should 
carefully take into account all the relevant dependen- 
cies of the perturbation integrals on the parameters of 
the theory. For this reason here the relevant perturba- 
tion contributions to the Hamiltonian vertex parameters 
are presented in a general form. In certain cases, the 
perturbation integrals are calculated for zero values of 
the external wave vectors, or, for zero values of certain 
Hamiltonian parameters [J [1] . In such cases some of the 
integrals give equal contributions to the respective RG 
equations, but in other cases, for example, the investi- 
gation of the RG stability matrix, one should carefully 
take into account all the differences in the contributions 
of the same integrals. The calculation of elements of the 
linearized RG stability matrix is made through the exact 
dependence of the perturbation terms on the parameters 
of the Hamiltonian and for this aim one needs to know 
the perturbation integrals in their initial general form. 
In all other aspects, the present consideration follows the 
standard prescriptions, described in Ref. [8,] and applied 
to the study of complex systems by means of RG and e- 
expansions in non- integer power of e (see, e.g., Ref. Q). 
In spite of the specific features of the present analysis, 
it does not contradict to the usual concepts and inter- 
pretations of e-series (see, e.g., Ref. [9[). In this paper 
some integrals appear for the first time in the RG the- 



ory of complex systems. These integrals together with 
known ones are evaluated to the accuracy needed for the 
RG analysis in two- loop approximation (See Sec. III). In 
RG studies in higher orders of the loop expansion some 
additional relevant terms of the same integrals should 
be calculated and used. In order to facilitate the repro- 
ducibility of the results and further investigations, some 
details of the calculations and an extended discussion of 
the most important stages of the RG investigation are 
presented (Sec. HI and IV). 

These remarks are important throughout the RG anal- 
ysis: (1) the calculation of fixed point (FP) coordinates, 
(2) the calculation of elements of the RG stability ma- 
trix, and (3) the calculation of critical exponents, includ- 
ing the stability ones. Besides, the derivation of the RG 
equations should be made with a considerable attention 
for reasons explained in Sec. HI. Note, that the RG inves- 
tigation can be performed in an alternative way, namely, 
the RG can be applied to a new effective Hamiltonian, 
which is a functional of the ■(/'-field only. This point is 
briefly discussed in Sec. II. 

The consideration of quantum effects exhibits for the 
first time an example of a complex quantum criticality 
characterized by a double-rate quantum critical dynam- 
ics (Sec. III.E). In the quantum limit (T — 0) the fields 
M and tp have different dynamical exponents, zm and 
z^, and this leads to two different upper critical dimen- 
sions: dff — 6 — z^ and di/ = 6 — zm- For this reason the 
more complete theoretical description of quantum effects 
on the properties of the zero-temperature (multi)critical 
point meets difficulties (Sec. III.E and Sec. V). The treat- 
ment of spin-triplet ferromagnetic superconductors with 
magnetic anisotropics {m < 3) and/or crystal symmetry 
lower than the cubic one requires a somewhat different 
RG analysis. These systems are considered in Sec. IV. In 
Sec. V the main results are summarized and discussed. 



II. EFFECTIVE HAMILTONIAN 

Relevant part of the fluctuation Hamiltonian of un- 
conventional ferromagnetic superconductors 0, 0] can be 
written in the form 



n 



E 

k 



(r + k^) \ij{k)\' + -{t + k^) 



\M{k)\' 



J2 M(fci).[^(fc2)x V*(fcl+fe2)] 
ki , Ala 



(1) 



where V ^ L is the volume of the d— dimensional sys- 
tem, the length unit is chosen so that the wave vector k 
is confined below unity (0<fc=|fc|<l),(7>0isa cou- 
pling constant, describing the effect of the scalar product 
of M and the vector product {ip x ip*) for symmetry in- 



dices n/2 ~ m = 3, and the parameters r = as{T — Tg) 
and t = af{T — Tf) are expressed by the critical tempera- 
tures of the generic {g = 0) superconducting (T,) and fer- 
romagnetic {Tf) transitions (as usual, the parameters as 
and Uf are positive). As mean field studies indicate 0,[6|, 
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Ts{P) is much lower than Tc{T) and Tf{P) ^ T/(P). 
As shown below the Hamiltonian (1) describes the main 
fluctuation efl'ects in a close vicinity of critical points in 
spin-triplet ferromagnetic superconductors. It is conve- 
nient to choose units, in which fc^ = 1, and the upper 
cutoff for the magnitude k = \k\ of the wave vectors 
in Eq. (1) is equal to unity. Some perturbation expan- 
sions within the RG investigation, in particular, those for 
isotropic systems (n/2 = m = 3) can be performed with 
the help of known representation £ja/3 of the components 
(-0 X of the respective vector product in Eq. (1). 

The fourth order terms (MMV'I'*, M^IVj^) in the total 
free energy (Hamiltonian) [y| have not been included 
in Eq. (1) as they are irrelevant to the present investi- 
gation. The simple dimensional analysis shows that the 
f/— term in Eq. (1) corresponds to a scaling factor &3-<i/2 
and, hence, becomes relevant below the upper border- 
line dimension du = 6, whereas fourth order terms are 
scaled by a factor as in the usual t/i'*— theory and 
are relevant below d < 4 (6 > 1 is a scaling number) 
Therefore, we should perform the RG investigation in 
spatial dimensions d — 6 — e, where the g-term in Eq. (1) 
describes the only relevant fluctuation interaction. Be- 
sides, the total fluctuation Hamiltonian 0] contains 
off-diagonal terms of the form kikjipa'tpp] i j and/or 
a p. Using a convenient loop expansion these terms 
can be completely integrated out from the partition func- 
tion to show that they modify the parameters (r, t, g) of 
the theory but they do not affect the structure of the 
model (1). So, such terms change auxiliary quantities, 
for example, the coordinates of the RG flxed points (FPs) 
but they do not affect the main RG results for the sta- 
bility of the FPs and the values of the critical exponents. 
Here these off-diagonal gradient terms will be ignored. 

The mean-field investigation of the total Hamiltonian 
shows that the f;-term in Eq. (1) triggers the spin- 
triplet superconductivity (M-trigger mechanism Q). 
The fourth order terms, mentioned above, do not play 
such a crucial role. They merely stabilize the ordered 
phases, provide a correct shape of the (P, T) phase do- 
mains for the specific crystal and magnetic symmetries. 
Neglecting these terms one looses the opportunity to 
describe the equilibrium order (M) = M — 6M and 
(■0) = ?/) — Sip, but as shown above, the main fluctuation 
effects in a close vicinity of the phase transition points 
can be totally taken into account. Here the attention 
is focussed on the behavior above the phase transition 
points where (M) — — 0, and hence, M = 6M and 
tp = 5ip. 

The simplest theory with upper critical dimension 
(i;7 = 6 is that of a classic scalar field Lp with 1^9^- 
interaction (see, e.g., Ref. [i|). Quite more complex 
problem, where dc/ = 6 and an expansions in powers 
of e = (6 — d) are used arises in the theory of tricritical 
Lifshitz points P, • In all these cases the upper criti- 
cal dimension djj is determined by a simple dimensional 
analysis in the so-called tree approximation [1] and/or 
by a check of singularities of the relevant perturbation 



integrals (See Sec. HI). The two mentioned examples are 
from theories of a single (vector or scalar) field whereas 
the present theory (1) describes two fields. Here the crit- 
ical dimension d^/ = 6 is a result of the simple fact that 
the total power of fields in interaction [g-) term in Eq. (1) 
is equal to three. Thus one may suppose that some of the 
results in this paper could be applied to critical phenom- 
ena in improper ferroelectrics |ll| . where interactions of 
two fields {(pi and tp2) of type ^pi^P2 ^-^so occur and the 
upper critical dimension is d[/ = 6. The present investi- 
gation will demonstrate, however, that for such interac- 
tion of two fields and upper critical dimension d[/ = 6, 
the RG analysis can lead to different results depending 
on the specific symmetry of the interaction term. The 
outcomes are two: a lack of stable FPs and physical ar- 
guments leading to a prediction of first order phase tran- 
sitions, or, the presence of a stable FP and, hence, a sta- 
ble critical behavior. While the particular symmetry of 
interaction term provides a new stable FP and, hence, a 
new critical behavior for symmetry indices n/2 = to = 3, 
one cannot be certain that the same prediction can be 
made for ferroelectrics without a specific investigation. 

One may consider several cases: (i) isotropic systems, 
namely, cubic crystal symmetry and isotropic magnetic 
order, when all field components ipa and Mj can be dif- 
ferent zero (n/2 = TO = 3), (ii) anisotropic systems when 
the total number (to + n/2) of field components is less 
than six. Note that the case (i) is of major interest to 
real systems, where fiuctuations of all field components 
are possible, despite the presence of crystal and mag- 
netic anisotropy that nullifies some of the equilibrium 
field components. 

The functional integral over the fields Mj{k) in the 
partition function Z — J DtpMexp{—H/T) can be ex- 
actly calculated, for example, by the total summation of 
perturbation series in powers of the interaction parameter 
g. In result, one will obtain a new effective Hamiltonian 
H{ip), which is a functional of a single field - 0. In this 
new V'-theory the effects of the magnetic (M-) subsys- 
tem will be "hidden" in the vertex parameters of H{ip). 
The effective Hamiltonian can also be treated by RG but 
this task is beyond the scope of the present consideration. 
The Hamiltonian (1) explicitly describes the fluctuations 
of the magnetization M and this is the advantage with 
respect to H{'ip). 

III. ISOTROPIC SYSTEMS 

A. RG equations 

Following Ref. [J, d, Q here we derive the Wilson- 
Fisher RG equations for isotropic systems (n/2 = to = 
3), described by the Hamiltonian (1), up to the two- 
loop order. The main results from the RG analysis of 
anisotropic systems {n/2+m < 6) can be obtained within 
the one-loop order and this point will be discussed in 
Sec. IV. The one loop contributions to the RG equations 
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FIG. 2: One loop diagrammatic contributions to the re- 
normalized correlation functions (ItAal^) and {|M-, |^), and to 
the g-term in Eq. (1); (•) represents the g-interaction, the thin 
legs represent Mj (fe), the tick legs with incoming and outcom- 
ing arrows stand for the field components tpaik) and V'a('s), 
respectively; the internal thin and thick lines represent the 
bare correlation functions { 1^/^0 1^)0 and {\Mjf)o, respectively. 



FIG. 3: Two loop diagrammatic contributions to the correlsr 
tion function {|i/)a(fe)|^). The symbols have been explained in 
Fig. 2. The arrows of the thick lines have been omitted. 




FIG. 4: Two loop diagrammatic contributions to the correla- 
tion function (|Mj(fc)|^). The symbols have been explained 
in Fig. 2. The arrows of the thick lines have been omitted. 




FIG. 5; A sum of two loop diagrams of order with a zero 
contribution to the renormalization of the g-term in Eq. (1). 
The symbols have been explained in Fig. 2. 




FIG. 6: Two loop diagrammatic contributions to the g-term 
in Eq. (1). The arrows of the thick lines have been omitted. 



are shown by diagrams in Fig. 2. The two loop dia- 
grammatic contributions to the renormaUzed correlation 

fimctions (|?A„(fc)p) and (|Afj(fc)|2) arc shown in Fie:. 3 
and 4, respectively. Although the perturbation theory of 
the Hamiltonian (1) is developed in a standard way, i.e., 
by an expansion in powers of the interaction parameter 
g, the derivation of the two-loop terms in the RG equa- 
tions is quite nontrivial because of the special symmetry 
properties of the interaction g-tcrm. For example, some 
diagrams with opposite arrows of internal lines, as the 
couple shown in Fig. 5, have opposite signs and compen- 
sate each other. The terms bringing contributions to the 
^-vertex are shown diagrammatically in Fig. 6. 

In Figs. (2) - (6) the thin external legs and thin internal 
lines, corresponding to the field components Mj{k) and 
the bare correlation function (|M,(A;)p)o, respectively, 
can always be supplied with incoming and/or outcoming 
arrows as this has been made for the thick legs and lines, 
representing the fields Va(fe) and V'aC')' ^^'^ hare 
correlation function {\ipa{k)\'^)o, respectively. The thin 
lines in Figs. (2) - (6) can take any orientation, because 
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the magnetization M{x) is a real vector and, hence, the 
Fourier amphtudes of the field components, Mj{k), obey 
the relation M*{k) = Mj{~k). In Figs. (2) - (6) the 
arrows of thin lines and legs are omitted but one should 



have in mind that in practical calculations with the help 
of these diagrams, arrows of any orientation can be used. 

The RG equations have been derived in the following 
general form: 



r' = b^-w {r - 2Jsf{r, t; 0)g^ - 2 [2Bi(r, t, 0) + 2B2(r, t, 0) + Bair, t, 0)] g^} (2) 
t' = b^-^- {t - 2 J,,(r, r; O)^' - 2 [A,{r, t, 0) + 4^2 (r, t, 0)] , (3) 

g' = ^3-d/2-„^-„M/2 1^ ^ ^)^3 ^ [^^(^^ ^ 4^2^^^ ^) ^ 2Z)3(r, t) + 2D4r, t)] g^} , (4) 

I 

6"* = 1 - 2a«/52 _ (25, + 262 + 63)5^ (5) Jss{r,r,k) = Jsf{r,r,k), (8) 



and 



fe-*" = 1 - 2assg^ - 2(ai + 202)/. (6) 

Here & > 1 is a scaling number, 77^, and rjM are the Fisher , ^ ^ 

exponents (anomalous dimensions of the fields tp and M, Jsir, t) = ^ — 2 2 — 2 ' 

respectively), and the perturbation integrals are given by J C^"^) (P + (P 

J,/ (r, fe) = I (p2 + [(p + fe)2 + ^] ' (7) 



A, (r- t fc) = / ^ ^^^ ^^ ^ (10) 

^^"^'"'"^ 7 (277)2^^ (p2 + ^)(p2+^)[(p^_p^)2 + i][(p^+fe)2+^][(p^ + fe)2+^]' ^ > 

(27r)2<i (pf + r)2 [(p^ + p,)2 + r-] [(pi + fe)2 + r] [(pi + 1] ' 



A2{r,t,k) I i^_^2d r«2 , ^\2 r/^ , „ \2 , ^1 r/^ , i,\2 , ^1 r/^? , ^1' (l^) 



Of fiA-[ 1 . . 

U^' ' ^ 7 (27r)2<^ (pf + r)2 [(pi + p2)2 + r] I(pi + fe)2 + [(p2 + ' I ^ 

Of f lA- (' d^Pid<'p2 1 , . 

'^"■'"'"^ J {2^r {pl+tn{p^+P2r+r][{p,+kr+r][{pl + ry 

/■' d'^pid'^p2 1 

3(^'*' )-y (27r)2^ (pf + r)(pi + t) [(pi - P2)2 + r] [{p, + kf + t] [(p^ + fe)2 + r] ' ^ ^ 



^ J {2nr<i (pi+t)ipi+t)ipi+r)mp,+P2r+rf' 



D2{r,t)= . „ , . . „ , . . „ , (16) 



(27r)2'^ (p2 + t){pl + t)(p2 + r)3 [(p, + p2)2 + r] ' 
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(27r)2rf [pI + tfipl + tfipl + r) [{p, + r] 



(17) 



Di{r, t) 



(18) 



r 



In Eqs. (5) and (6), a^/, Oss, ai, 02, 61, 62 and 63 are dif- 
ferences of integrals, given by the rule 



(19) 



fc=0 



where Yi denotes one of the integrals Jsf , Jss , Aj or , 
and yi (r, i) stands for = a/ or 5/ and accordingly, the 
index I denotes either the symbols s/, ss or the numbers 
j = 1, . . . ,3. In Eqs. (7)-(18) / denotes an integration 
in restricted domains of wave numbers [l|, for example, 
for the integral Ai this domain is given, as implied by 
the terms in the denominator of the integrand, by the 
inequalities 1/6 < _Pi < 1, l/b < \pi — P2I — 1/^ — 
|Pi + fe| < Ij where the external wave vector fc obeys the 
condition < fc < l/b, and i = 1,2. 

The general RG Eqs. (2) - (6) are the starting point 
of the present RG analysis in two loop approximation. 
The same Eqs. (2) - (6) can be considered as reliable 
stage in the derivation and analysis of RG equations in 
higher orders of the loop expansion. The integrals (7) - 
(18) correspond to the most general diagrammatic rep- 
resentation of the respective perturbation terms given in 
Figs. 2-6 with the only difference that the dependence 
of the integrals Dj{r, t, ki, ^2) on the external wave vec- 
tors ki and ^2 has been neglected as irrelevant for the 
RG treatment [l|, namely, Dj{r,t) in Eqs. (15) - (18) 
denotes Dj{r,t, 0,0) for any j — 1,...,4. For d ~ 6 



the integrals J3(0,0), Dj {0,0) exhibit infrared logarith- 
mic divergencies at 6 = cxd, which means that the upper 
critical dimension is djj = 6 and the e-expansion should 
be developed in powers of e = {6 — d). This conclusion 
is in a total conformity with the dimensional analysis, 
mentioned in Sec. II. 

In the framework of the two loop approximation cer- 
tain simplifications of Eqs. (2)-(6) can be made with- 
out any effect on the final results of the RG analysis 
[E H, g. One can set Dj{r,t) « Dj{0,0) = Dj and, 
hence, D2 = -D3 ^ Ci{h), and Di = D/^ = C2{b), where 



and 



Ci{b) 



C2{b) 



' d'^pid'^P2 



(27r)2rf plplip,+P2)' 



' d-^pid-^P2 



(271)2-^ pfp|(Pi+P2)'' 



(20) 



(21) 



Moreover, setting r = i = in Eqs. (5) and (6) one 
obtains that 

(261-^262-^63) = (4a2 + ai). (22) 
For the needs of the two loop RG analysis, the relevant 
6-dependence of the integrals (7) - (18) has been obtained 
in the following form: 



J,f = [(2 + e)(l - 6-2) - 4(r + t)\nb - 2e6-2ln6] 



(23) 



Kd-i 
32 



2 + 1^ ) ln6 - (b' ~ l)(r + i) + e(ln6)' 

TT 



(24) 



where 



/ dy sin^'y [In (siny)] (4cos^2/ — l) , (25) 
Jo 

ass(r, t) = asf{r, r), (26) 



ai = -^^N + 2(ln6)2] (28) 



K3 



a2 = 



3072 



i [Illn6 + 6(ln6)2 -9(6^ - 1)] , (29) 



J3 = ^ [21n6 -t- e(ln6)2 - (6^ - l)(2r + t)] (27) 



Ci = W - 1) - 111^^' - 6(ln6)^] , (30) 
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C2 = 5^|^[ln& + 2(ln6)2], (31) 

where Kd ^ 2^~'^tt-'^/'^ /T{d/2). Note the following useful 
formulae for r = t = 0: (0,0,0) = -83(0,0,0), where 

Ai{0,0,0) = -^KdKd-i, (32) 
(0,0,0) = Bi (0,0,0) ^2(0, 0,0), where 



^2(0,0,0) = —i^dXrf_iln6, (33) 



^ = -4C2, ^^-3C,-C„ (34) 
or or 



dA, 3M._,^^^^ '-^^-C,, (35) 



dt 



16 ' dt 



^^-2C2-Ci, (36) 



dr 



dr 



f)H 3 
— ^ - -2C2 - —KaKd-i\nb, 
Or 16 



— — Ci — C2 



and 



^ = -2C2. 



dt 



dt 



(37) 



(38) 



(39) 



B. One loop approximation 

Neglecting the g'^-terms in Eqs. (2) and (3) as well as 
the g^-terms in Eq. (6), using the respective integrals Usf 
and Sss as given by Eqs. (24)- (26) to the zero approxi- 
mation in e = (6 — d), r, and t, we obtain that 



Kd-l *2 



Using the expansion 6^ « 1 -I- ln6, and bearing in mind 
that Ke = l/GAn^ and K5 = l6Ke/3, i.e., that K5 = 
16Kq/3, one easily obtains that Eq. (42) has the simple 
form 



eg* = 0. 



(43) 



(40) 



For e 7^ this equation yields only a Gaussian FP (GFP) 
{g* = go — 0, re ^ to ^ 0), which is stable for di- 
mensions d > 6 but is unstable for d < 6. As usual the 
GFP describes exponents, corresponding to the Gaussian 
model for d < 6 and mean-field exponents for d > 6. For 
example, one may use Eqs. (2)-(4) to obtain that the 
correlation length exponents i'^ and i>m, corresponding 
to the correlation lengths of the ip- and iW-subsystems, 
respectively, have no e-corrections {v^ = vm — 1/2) i f^nd 
that the stability exponent, corresponding to the param- 
eter 5, is given by j/g = e/2. 

Further, Eq. (43) shows the first particular property 
of the present RG analysis which has no analog in other 
systems. For d = djj = Q, Eq. (43) allows an arbitrary 
value of g* within the framework of usual physical con- 
ditions. This means that for d = 6 a new FP exists and 
is characterized by a real coordinate g* > 0, which has 
no specified value and can take any positive real value. 
This object can be called arbitrary FP (in short, AFP). 
Note, that complex and negative real values of g* are 
outside the physical domain of values of parameter g. 
FPs with such coordinates are usually called unphysical 
(UFP) [9]. The exponents ry,/, and rim and the FP co- 
ordinates rA and tA are given by Eqs. (40) and (41), 
where g* € (0,oo). The exponents corresponding to 
AFP can be calculated in a standard way mi. The re- 
sults are: 1^1 = (l/2-53V3847r3), z/2 = {1/2 + /96n^), 
and yg — 0. It is seen that in one loop approximation 
AFP has a marginal stability {yg = 0). This means that 
the stability properties of this FP should be investigated 
in two loop order of the theory but this is beyond the 
aims of this consideration. AFP exists only at dimen- 
sions d = 6 and, hence, it has no real physical signifi- 
cance. Rather, such FPs are of pure academic interest. 
Thus, there is no real motivation for its further investiga- 
tion. In the remainder of this paper only FP with some 
physical significance and stability in dimensions d < 6 
will be considered. 



where g* is any FP value of the vertex parameter g. The 
one loop FP values r* and t* are obtained by Eqs. (2) and 
(3) where one should neglect the ^'^-terms. The result is 



Kdg^ 



(41) 



Finally, neglecting the g'^-terms in Eq. (4) the equation 
for the possible FP values g* of g is obtained. Taking 
the integral J3 from Eq. (27) to zeroth order in e, r, and 
t one obtains that g* is determined by the zeros of the 
following simple equation: 
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^(e-3r,)/2 ^_ g*3Kdlllb) 



(42) 



C. Nontrivial FP in two loop approximation 

With the help of Eq. (22), (24) -(26), (28) and (29), 
Eqs. (5) and (6) can be solved with respect to the expo- 
nents 77^ and rjM- One obtains that the respective equa- 
tions for these exponents are identical and, hence, the 
solutions should be equal 77 = 77^ = rjM- In the respec- 
tive equation for the exponent 77, the terms containing 
factors (6^ — 1), see Eqs. (24) and (29), compensate one 
another and only terms containing 6-dependent factors of 
type ln6 and (ln6)^ remain. Expanding 77 to fourth order 
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in g, namely using 



T] = aig + a2g 



*4 



(44) 



the following equation for the coefficients ai and a2 are 
obtained 



-»? + 



1 2 , KeKd-i 2 

-g + g - 



768 



lnb + 



eKd 



1 2 Ctl 4 

q q 

16 ^ 2 ^ 



-1 „4 



128 



/ (ln6)^ = 



(45) 



r 



(hereafter in this Section the superscript (*) of g* will be 
often omitted). The Eq. (45) is an expansion in powers 
oi g^ and contains terms of types: g^\nb, g*lnb, g'*(ln6)^, 
eg^lnb, and eg'^{\nb)'^. These terms can be grouped in 
three sums of terms of three types: .g^lnb, g^lnfe, and 
5^(ln6)^. Having in mind that Eq. (45) should be satisfied 
for any 6 > 1 , all sums of the mentioned types should be 
equal to zero. Taking into account Eq. (44) and setting 
the sum of terms of type gHnb in Eq. (45) equal to zero, 
one easily checks the one loop result ai = Kd-ig*'^ /8; c.f. 
Eq. (40). Now one should set the sum of terms of type 
g^liib equal to zero. This yields the following equation 
for g: 



a2 + 



768 



g = 



KKd 



(46) 



Now putting the sum of terms of type (ln6)^ equal to zero 
one finds 



64 



Kd-i 



-€. 



(47) 



Eq. (47) reproduces the one loop result for ai, only if the 
r.h.s term is small compared to the l.h.s. ones, namely, 
if g^ ~ e*', where y < 1. The shape of Eq. (46) is in a 
conformity with this assumption. In fact, both Eqs. (46) 
and (47) show that the opposite supposition (y > 1) does 
not allow the existence of FPs of type g* > 0. Accept- 
ing the hypothesis y < I, and neglecting the r.h.s term 
in Eq. (46), one gets 0:2 = — 13i<'J_i/768, and, hence, 
Eq. (44) takes the form 



Kd-i 2 
V= —5—9 - 



768 



-^g\ 



(48) 



where g = g* > 0. 

In order to find the FP values 17* ~ e** > of parameter 
g, the RG Eq. (4) should be investigated. Using the 
relations between the integrals £>j(0,0) and Ci and C2, 
as discussed in Sec. III.B, this equation takes the form 
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(49) 



Setting in Eq. (49) g' = g = g* and the values of integrals 
J3, Ci and C2, given by Eqs. (27), (30), and (32), respec- 
tively, one obtains an equation for g > 0, which contains 
terms with three different 6-dependent factors: {b^ — 1), 



ln6 and (ln&)^. The two terms with factors (6^ — 1) cancel 
each other, the sum of terms with (ln6)-factors become 
equal to zero, as should be, if the nonzero FP value of g 

is given by 



8(37r3)'/'(2e/13) 



1/4 



(50) 



The terms of type g'^{lnb)'^ compensate one another, pro- 
vided Eq. (50) for the FP value of g takes place. Note, 
that the assumption g"^ ~ with y < 1 leads to eg^ <C g^ 
and, therefore, terms of order eg^ ^ gi+y and order in 
the FP equation for g are small and can be safely ignored. 

Thus the assumption g*'^ ^ yields an essentially new 
nontrivial FP (50) for y = 1/2. In the framework of this 
new scheme of RG investigation the two loop order gives 
RG results up to the first order in e whereas the one loop 
order can be used for calculation within an accuracy of 
order e^/^. 

With the help of Eqs. (2) and (3) as well as Eqs. (8), 
(23), (32), (33) and the relation between the integrals Bj 
and Ai and A2, one can easily obtain the FP coordinates 



Kdg^ + 



g 



487r6 

r* =t*, or, using Eq. (50), 



+ 0{g%-Hnb), 



r* =t* = 3\ 



'2e 1536 
13 + ^" 



(51) 



(52) 



D. Critical exponents in two loop approximation 

Using Eqs. (48)and (50) the critical exponent rj can be 
written in the form 




(53) 



The critical exponents and um of the correlation 
lengths of magnetic and superconduction subsystems, re- 
spectively, as well as the stability exponent yg describing 
the stability of the FP (50) with respect to the inter- 
action parameter g, can be obtained as eigenvalues of 
the linearized stability matrix /t = (9/u-/9/ij) of the RG 
transformation (2) - (4); here /x = ^^2, ^^^) = {f'l't^g) 
is a notation of a vector in the parameter space (r, t, g) 
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of the Hamiltonian (1). Following Refs. [U, H, [1], the 

matrix elements — (dr'/dr), /ii2 = {dr' /dt), = 
(dr' /dg),. . . , are obtained in the form: 

(54) 

where x = Kj^g^, 



112 = 2.t62-''+''^/3 <^ Infc 



(62 _ 1) _ ln6 



Ml3 — ^23 
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;i2i = 4x62-''+^ J Infe 



M22 



/2-j)+4a:-^ 



3(&2 - 1) - -ln6 



[1 + 4x2(lnfe)2 



Mai = 



1-62 



32 



(55) 



(56) 



(57) 



(58) 



(59) 



and 



A3 = &^-^-^(ln6)2e, 



where the exponents 



2/2 = 2- 



197 



and 



2/3 



(63) 



(64) 



(65) 



(66) 



have been identified in the following way: yi = yr, 
y2 = yt and y^ = yg. The negative sign of yg for c? < 6 
means that the nontrivial FP (50) is stable and describes 
a critical behavior. 

The correlation length critical exponents = l/yr 
and vm — ^/yti corresponding to the fields and M, 



and 



M33 



^ 5(e-3i))/2+a:-91a:V' 



36 



(62-l)-2(ln6)2 



Ail 



(60) 
is 



The solution of the eigenvalue equation |/t 
quite simple for GFP, where = x = 0. But the solution 
of the same problem for the nontrivial FP, given by Eqs. 
(50) and (52), is obtained by quite lengthy calculation. 
Here the main steps of this calculation will be outlined. It 
seems convenient to emphasize that this treatment does 
not depart from the calculational schemes known from 
preceding papers fl, l9i| . 

Having in mind that x = Kdg^ and -q in Eqs. (54)-(60) 
can be represented by e as given by Eq. (50) and (53), 
the coefficients of the algebraic equation of third order 
for A are calculated as functions of e. The three roots 
(Ai,A2,A3) of the eigenvalue equation are calculated in 
powers of e^/^; 



A = /3o + /3ie'/' + /32e 



(61) 



The calculation of the eigenvalues Aj = Aj(6)6^^ of the 
matrix jl should be performed very carefully, in particu- 
lar, by focussing a special attention to the behavior of the 
large dangerous terms of type 6^ and 62(ln6), (6 1) [1| 
in the elements of the matrix /x. In the course of 
the calculation the most part of these terms cancel one 
another and, hence, have no effect on the final results for 
the critical exponents. As usual Q, some of these terms 
produce large scaling amplitudes Aj{h). As in the usual 
(/>*-theory the amplitudes Aj depend on the scaling 
factor 6. The result for the eigenvalues is 



Ai,2 = 6^1'^ 1 



^62, 
13 



(62) 



1 „ /2e 5e 
VM = - + 2\ . 

2 V 13 156 



(67) 



(68) 



These exponents describe quite particular multi- 
critical behavior, which differs from the numerous exam- 
ples known so far. For d = 3, = 0.78, which is some- 
what above the usual value, v ^ 0.6-^0.7 near a standard 
phase transition of second order [l[ but fM = 1.76 at the 
same dimension (d = 3) is unusually large. The fact 
that the Fisher's exponent [l| r; is negative for d = 3 
does not create troubles because such cases are known 
in complex systems, for example, in conventional super- 
conductors |12|. Perhaps, a direct extrapolation of the 
results from the present e-series is not completely reli- 
able because of the fact that the series has been derived 
under the assumptions of e ^ 1 and under the conditions 
e^/2& < 1^ ei/2(ijj5) ^ X, provided 6 > 1. These condi- 
tions are stronger than those in the usual ^''-theory P, . 
Using the known relation fl| 7 = {2 — ifjv, the suscepti- 
bility exponents for d = 3 take the values 7^ = 2.06 and 
7m = 4.65. These values exceed even those correspond- 
ing to the Hartree approximation [ll| (7 = 2i/ = 2 for 
d = 3) and can be easily distinguished in experiments. 

Note, that here the interpretation of the e-series and, in 
particular, extrapolations of e-results up to e = 3 should 
be made with a caution and according to the remarks 
presented in Ref. Q. The extrapolations of results for 
critical exponents, for exainple, Eqs. (67) and (68), to fi- 
nite values of e are reliable [9| , if and only if e-terms give a 
small correction to the value of the respective exponent. 
For example, the e-corrections in Eq. (68) do not satisfy 
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this rule for e = 3 and, hence, the result for vm for d = 3, 
given by Eq. (68), seems unreliable. But the Eq. (68) is 
an useful result, which indicates that the value of the 
exponent -qj^i in real three dimensional systems, is large 
compared to that given by usual theories and this is a 
reliable conclusion. Let us emphasize that this relatively 
large value of the exponent vm has a physical explana- 
tion in the fact the g-ievva includes the first-order power 
of M . The relatively low power of M in the (/-term in 
Eq. (1) is the reason for unusually strong magnetic fluc- 
tuation effects on the critical behavior. These effects are 
responsible for the relatively large values of the critical 
exponents, in particular, for the large values of the ex- 
ponents VM and 7Af. This example has been given to 
show the way, in which the present results should be in- 
terpreted in discussion of real three dimensional systems. 



E. Notes about the quantum criticality at zero 
temperature 

The critical behavior discussed so far may occur in a 
close vicinity of finite temperature multi-critical points 
(Tc — Tf>{)) in systems possessing the symmetry of the 
model (1). In certain systems, as shown in Fig. 1, this 
multi-critical points may occur at T = 0. In the quantum 
limit (T — > 0), or, more generally, in the low-temperature 
limit [T ^ /i; /I = (i, r);kB = 1] the thermal wavelengths 
of the fields M and ip exceed the inter-particle interaction 
radius and the quantum fluctuations become essential for 
the critical behavior 0, [l^ . 

The quantum effects can be considered by RG analysis 
of comprehensively generalized version of the model (1), 
namely, the action S of the referent quantum system. 
The generalized action is constructed with the help of the 
substitution {-H/T) S[M{q),i}{q)]. The description 
will be given in terms of the (Bose) quantum fields M{q) 
and Tpiq), which depend on the {d + l)-dimensional vec- 
tor q ~ (ujijk); LUi = 2-kIT is the Matsubara frequency 
(/i = 1;Z = 0,±1,...). The fc-sums in Eq. (1) should 
be substituted by respective g-sums and the inverse bare 
correlation functions [r + k"^) and {t + k"^) in Eq. (1) con- 
tain additional w/— dependent terms, for example [3. ITsj 

{\M<l)?r' ^\^i\ + k^+r. (69) 

The bare correlation function {\Mj{q)\'^) contains a term 
of type \LOi\/k^ , where 9 = 1 and = 2 for clean and dirty 
itinerant ferromagnets, respectively [isj . 

Quantum dynamics of the field '0 is described by the 
bare value z = 2 of dynamical critical exponent z = 
whereas the quantum dynamics of the magnetization cor- 
responds to zm = 3 (for = 1), or, to zm = 4 (for 
9 = 2). This means that the classical-to-quantum di- 
mensional crossover at T ^ is given by d — > (d -\- 2) 
and, hence, the system exhibits a simple mean field be- 
havior for d> A. Just below the upper quantum critical 
dimension = 4 the relevant quantum effects at T = 



are represented by the field ip whereas the quantum {uji- 
) fluctuations of magnetization are relevant for d < 3 
(clean systems), or, even for d < 2 (dirty limit) p^ . 
This picture is confirmed by the analysis of singularities 
of relevant perturbation integrals. Therefore, the quan- 
tum fluctuations of the fleld ij: have a dominating role for 
spatial dimensions c? < 4. 

Taking into account the quantum fluctuations of the 
field ijj and completely neglecting the ^/-dependence of 
the magnetization M , cq = (4 — (i)-analysis of the gen- 
eralized action S has been performed within the one- 
loop approximation (order ej). In the classical limit 
{r/T ^ 1) one re-derives the results, already reported 
above together with an essentially new result, namely, 
the value of the dynamical exponent = 2— (2e/13)^/^, 
which describes the quantum dynamics of the field il). In 
the quantum limit (r/T 3> 1, T — > 0), static phase transi- 
tion properties arc affected by the quantum fluctuations, 
in particular, in isotropic systems {n/2 = m = 3). In 
this case, the one-loop RG equations, corresponding to 
T = 0, are not degenerate and give definite results. The 
RG equation for g, 

9' = b^"/^g (l + ^Info) , (70) 

yields two FPs: (a) a Gaussian FP {go = 0), which is 
unstable for d < 4, and (&) a FP {g'^)* = -Un^eo, which 
is unphysical [{g^)* < 0] for d < A and unstable for d > 4. 
Thus the new stable critical behavior, corresponding to 
T > and c? < 6, disappears in the quantum limit T 0. 

At absolute zero (T = 0) and any dimension d > the 
P— driven phase transition (Fig. 1) is of first order. This 
can be explained as a mere result of the limit T 0. 
The only role of the quantum effects is the creation of 
the new unphysical FP (&). In fact, the referent classical 
system described by TC from Eq. (1) also looses its stable 
FP (8) in the zero-temperature (classical) limit T — s- 
but does not generate any new FP, because of the lack 
of g'^-teim in the equation for g'; see Eq. (13). At T = 
the classical system has purely mean field behavior [2], 
which is characterized by a Gaussian FP {g* = 0) and is 
unstable towards T-perturbations for < d < 6. This 
is a usual classical zero temperature behavior, where the 
quantum correlations are ignored. For the standard 0^— 
theory this picture holds true for d < 4. 

One may suppose that the quantum fluctuations of 
the field are not sufficient to ensure a stable quantum 
multi-critical behavior at Tc = Tp = and that the lack 
of such behavior is a result of neglecting the quantum 
fiuctuations of M. One may try to take into account 
these quantum fluctuations by quite special approaches 
from the theory of disordered systems, where additional 
expansion parameters are used to ensure the marginal- 
ity of the fluctuating modes at the same borderline di- 
mension djj. It may be conjectured that the techniques, 
known from the theory of disordered systems with ex- 
tended impurities, cannot be straightforwardly applied 
to the present problem and, perhaps, a completely new 



11 



approach should be introduced. 



IV. ANISOTROPIC SYSTEMS 

For anisotropic systems, mentioned in Sec. II as case 
(ii), the RG analysis leads to different results. The 
RG equations have no stable FP for dimensions d < 6 
and, therefore, one may conclude, quite reliably, that 
the Afj'0Q,'(/'^-interactions in the Hamiltonian (1) induce 
fluctuation-driven first order phase transitions. The RG 
result for a lack of FP is not enough to make conclu- 
sions about the order of the phase transition, but in our 
case there are additional strong heuristic arguments as 
well as arguments from preceding mean-field investiga- 
tions 0. The heuristic argument is that the account 
of anisotropy effects often leads to change of the phase 
transition order from second order phase transition, when 
such effects are ignored to a first order phase transition, 
when the anisotropy is properly taken into account (see, 
e.g., Ref. Q). Besides, a recent mean field study shows 
that first order phase transitions and multi-critical points 
are likely to occur in both isotropic and anisotropic fer- 
romagnetic superconductors for broad variations of the 
theory parameters {r^t,g, ...). These non-RG arguments 
strongly support the present point of view that the lack 
of FP of RG equations for anisotropic systems can re- 
liably be interpreted as indication for fluctuation-driven 
first order phase transitions in dimensions d < 6 and, in 
particular, for c? = 3. 

Now, one has to prove that RG equations, correspond- 
ing to anisotropic systems do not exhibit stable FP for 
d < 6. The case (ii) of anisotropic systems, mentioned 
in Sec. II, contains several subclasses of systems. For ex- 
ample, one may have cubic crystal anisotropy with XY 
magnetic symmetry (Mi,M2,0), i.e., magnetic symme- 
try index m = 2, or, magnetic anisotropy of Ising type: 
(0,0, M3), i.e., m — 1. Another example is the tetrago- 
nal crystal anisotropy, where the symmetry index for the 
superconducting order is n/2 = 2, i.e., {^pi,ip2, 0) and the 
magnetic symmetry can be of two types: XY order sym- 
metry (m = 2) and Ising symmetry (m =1). These and 
other examples of anisotropic systems with total symme- 
try index (n/2 -I- m) < 6 can be considered separately to 
prove the lack of stable FP in each of the cases. 

We shall demonstrate the lack of stable FP for d < 6 
for uniaxial (Ising) magnetic anisotropy (0, 0, A/3) and 
cubic crystal symmetry (n/2 = 3). In this case, the one 
loop RG equations have the form 

r' ^b'-^[r~Jsfir,t,0)], (71) 
t' ^b^-^'^[r-2J,s{r,r,0)], (72) 

and 



where the exponents rji (= ry), r]2 (= r]) and tjm corre- 
spond to the fields components V'l, V'2, and M3, respec- 
tively. The exponents ?/i and 772 are equal because the 
equations for them are the same 

b^'-- ^l-a^fg^. (74) 

From Eq. (74) one obtains 771 = 772 = ?y = K^-ig^ /16- 
Because of this equality of 771 and 772, Eq. (71), which 
describes the RG transformations of the two Hamiltonian 
terms of type r|7/;ip and r|i/;2p, is the same. This fact 
ensures a self-consistency of RG transformation for the 
parameter r. The one loop equation for the exponent 77M 
is given by Eq. (7), provided the g'^-teiui in the r.h.s. of 
this equation is ignored. Thus one finds that 77M = 2rj = 
Kd-ig^ as is in isotropic systems. 

It should be emphasized that the perturbation series 
does not give self energy contributions to the |?/;3p-term. 
That is why the RG transformation is performed after 
a simple procedure of integration out of this field in the 
functional integral for the system partition function. The 
integration is exact as the respective integral is Gaussian. 
The effect is an irrelevant contribution to the system free 
energy. Thus, the field component t/js is redundant in 
this consideration and does not participate in the RG 
transformation. On the other hand, if one keeps the field 
component ip in this RG transformation and considers it 
at the same footing as the other fields (t/'i, 7/12, M3) the 
resulting RG transformation of the term (r -I- /c^)|'03(fc)p 
in Eq. (1) will produce a transformation for the param- 
eter r, r' — b^r, which contradicts to the relation (71). 
Hopefully, this is not the right way of RG treatment and 
here its discussion should be understood only as a note 
of caution. 

Using Eq. (73), Eq. (27) for the integral J3 to zeroth 
order in e, r, and t, as well as the above results for the 
exponents rj and 77^, one can easily obtain FPs. The 
GFP {go — 0) exists and is stable for d > 6. The second 
FP is given by g*^ — — 967r^e, which means that this FP 
is unphysical for e > 0, i.e., for dimensions d < 6. For 
d > 6 this FP is physical, i.e., g* has a real positive value 
but for these high dimensions it is unstable towards the 
parameter g (jjg — — lle/2 > for e < 0). In this way we 
proved the lack of stable FP for dimensions d < 6. The 
result cannot be changed in the next orders of the loop 
expansion. 

The RG analysis, performed above, is valid also for 
systems with uniaxial magnetic anisotropy (0, 0, M3) and 
tetragonal crystal symmetry, (t/^i, -02, 0). The only differ- 
ence is that in these systems the redundant field compo- 
nent -03 is equal to zero and one does not need to perform 
a functional integration over redundant fields. 

For tetragonal symmetry and bi-axial (XY) magnetic 
anisotropy (Mi, M2, 0), the term in Eq. (1) is equal to 
zero and, hence, Eq. (1) describes a simple Gaussian fluc- 
tuations. In this case one must consider the fluctuation 
effects coming from the fourth order terms in the general 
effective Hamiltonian 0, Q . These terms may lead to the 
appearance of familiar types of FPs, which describe bi- 
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critical and tertracritical points of phase transitions (see, 
e.g., Ref. d). 

A reliable conclusion for the fluctuation-driven first or- 
der phase transition in systems with cubic crystal sym- 
metry and XY magnetic order {M3 = 0) can be made 
only after a check by the one loop RG investigation, but 
the results established so far imply that the appearance of 
stable FP in such systems is quite unlikely. The common 
difference in the results from one loop RG equations for 
isotropic and anisotropic unconventional ferromagnetic 
superconductors is in number factors and scaling expo- 
nents of type b'^ in the respective RG equations. The new 
stable FP (50) occurs as a result of quite special set of 
number coefficients and scaling factors in the RG equa- 
tions, which is not the case in anisotropic systems; c.f. 
the one loop order terms in Eqs. (2)- (6) on one hand and 
the terms in Eqs. (71)-(74) on the other hand. 

V. CONCLUSION 

The general RG equations for ferromagnetic supercon- 
ductors with spin-triplet Cooper pairing were derived and 
analyzed up to the second order in the loop expansion. 
For cubic crystals with isotropic magnetic order a new 
universality class of critical behavior was predicted. The 
main features of this new critical behavior were estab- 
lished and the critical exponents were calculated. It has 
been shown that lower crystal and magnetic symmetry 
produces different fluctuation effect - a fluctuation change 
of the phase transition order. The way of interpretation 
of the results and the extrapolation of e-series to real 
dimensions {d — 3) have been discussed at the end of 
Sec. III.D. The strong fluctuation interactions of type 
Mjipatpi3 have a crucial effect on the quantum criticality 
at zero temperature and some new features of this quan- 



tum phase transitions have been outlined in Sec. III. E. 
The complete understanding of this quantum phase tran- 
sition requires further theoretical investigations and, per- 
haps, some new ideas of calculation. As mentioned in 
Sec. III.E, the satisfactory consideration of the quantum 
fluctuations of both fields M and ijj requires a new RG 
approach, in which one should either consider the differ- 
ence (zm — z^) as an auxiliary small parameter or in- 
vent a completely new theoretical scheme of description. 
This problem is quite general and presents a challenge 
to the theory of quantum phase transitions 0. Within 
this research the present author has not been able to 
give a comprehensive solution of the problem and, hence, 
the discussion of quantum effects presented in Sec. III.E 
should be accepted as a preliminary outline of general 
problems, rather than a report of a complete description 
of this type of quite complex quantum phase transitions. 

The results can be of use in interpretations of recent 
experiments in UGe2 , where the magnetic order is uni- 
axial (Ising symmetry) and the experimental data, in ac- 
cord with the present consideration, indicate that the C- 
P phase transition is of first order. Systems with isotropic 
magnetic order are needed for an experimental test of the 
newly described multicritical behavior. The present re- 
sults can be applied to any natural system of the same 
class of symmetry, although this report is based on par- 
ticular example of itinerant ferromagnetic compounds. 
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